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The Schramm-Loewner evolution (SLE) describes the continuum limit of domain walls at phase 
transitions in two dimensional statistical systems. We consider here the SLE in Z{N) spin models 
at their self-dual critical point. For N = 2 and A'' = 3 these models correspond to the Ising and 
three-state Potts model. For A'^ > 4 the critical self-dual Z{N) spin models are described in the 
continuum limit by non-minimal conformal field theories with central charge c> 1. By studying the 
representations of the corresponding chiral algebra, we show that two particular operators satisfy a 
two level null vector condition which, for > 4, presents an additional term coming from the extra 
symmetry currents action. For N = 2,3 these operators correspond to the boundary conditions 
changing operators associated to the SLE16/3 (Ising model) and to the SLE24/5 and SLE10/3 (three- 
state Potts model). We suggest a definition of the interfaces within the Z{N) lattice models. The 
scaling limit of these interfaces is expected to be described at the self-dual critical point and for 
iV > 4 by the SLE4(jv+i)/(jv+2) and SLE4(jv+2)/(jv+i) processes. 



I. INTRODUCTION. 



The Schramm-Loewner evolutions (SLEs) are random growth processes that generate conformally invariant curves 
in two dimensions (2D). SLEs yield probability measures on the continuum limit of non-crossing interfaces in 2D 
statistical lattice models at criticality with conformal invariance. SLEs has been proved successful to a more complete, 
and in some case mathematically rigorous, description of fractal curves in critical percolation loop erased walk 
level lines of height models [sl] and domain boundaries at phase transitions (see e.g. and references therein). The 
SLEs consider directly the geometrical characterization of non-local objects in 2D critical systems and complement 
the powerful tools provided by the conformal field theories (CFTs) techniques. 

On the other hand, CFTs focus on correlation functions of local operators which are the scaling limit of lattice 
variables. A first family of CFTs, the minimal CFTs, calculates these correlation functions by studying the infinite 
constraints imposed by the conformal invariance in 2D systems. The Hilbert space of these theories is constructed 
from the highest weight representations of the Virasoro algebra formed by the generators of the conformal symmetry. 

The relation between SLEs and CFTs has been worked out in 0,0, 0|- In this respect, an important role is played 
by the boundary conditions changing operators (b.c.c.) which generate the boundary conditions from which the 
curve grows [1^ . [isj . The key property of these operators is to satisfy particular relations under the action of the 
conformal symmetry generators which lead to second order differential equations for their correlation functions. The 
link between the SLEs and CFTs is derived by comparing these second order equations to the Fokker-Plank equations 
coming from the Brownian motion driving the SLE. 

The SLEs/CFTs connection is well established in the case of minimal CFTs, which have a central charge c < 1. 
However, many critical models in condensed matter and statistical physics posses, in addition to the conformal 
invariance, symmetries in some internal degree of freedom, such as the SU (2) spin-rotational symmetry in critical 
quantum spin chains Q or the Zpf symmetry in spin lattice models jiol. fll|. In the continuum limit, these additional 
symmetries are enhanced by the presence of chiral currents which form, together with the energy-momentum tensor, 
more structured algebras which present the Virasoro one as a sub-algebra. The space of local fields occurring in a 
CFT with additional symmetry, the non-minimal CFT, corresponds to the representations of the associated chiral 
algebra. These theories have in general a central charge c > 1. 

The connection between SLEs and non- minimal CFTs has been first addressed in where the relation between 

stochastic evolutions and superconformal field theory was investigated. More recently, the connection between SLE 
and Wess-Zumino-Witten models, i.e. CFTs with additional Lie-group symmetries, has been studied by very different 
approaches [lBili3- particular, it has been shown in that a consistent SLE approach is possible provided that 
an additional stochastic motion in the internal symmetry group space in introduced. For an SU (2) symmetry, for 
instance, the SLE describes a process which carries a fluctuating additional spin 1/2 degree of freedom. 
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Inspired by these results, we discuss in this paper a possible SLE approach to CFTs with additional Zn symmetry. 
The conserved currents associated to the extra symmetry, called parafermions, have fractional spin. We consider the 
first of such parafermionic theories with central charge c = 2{N — l)/(A^ + 2) 10]. The Z{2) and Z{3) models coincide 
with some particular minimal CFT where the domain wall boundary conditions and the associated b.c.c. operator are 
known. We exploit these pieces of information to identify the possible b.c.c. operators for general Zn symmetry. We 
show that, under the action of the parafermionic conserved currents, these operators satisfy algebraic relations similar 
to the ones found in the case of minimal models. From these CFT predictions and in analogy with the approach 
proposed in p^j . we introduce an additional stochastic motion driven by the action of the parafermionic currents. 

The strong motivation to study these theories lies on the fact that they describe the continuum limit of critical 
lattice models, the iV-states generalization of the Ising model which reduce to the Ising and three state Potts model 
for N = 2 and = 3. This allows for the identification of the SLE interfaces on the lattice and opens the possibility 
of a direct numerical verification. 

The paper is organized as follows. We review the connection between SLE and minimal CFTs in Section |TT1 Section 
mil provides the extension to CFTs with extra Lie symmetries. In Section [Tvl after an introduction to the operator 
content of the parafermionic theories, we define a stochastic motion in the Zn internal space. Further, by studying the 
parafermionic algebra representations, we derive some algebraic relations related to the b.c.c. operators. In Section 
|V]we shall discuss the definition of interfaces on the lattice. We conclude the paper in Section IVll 



II. STOCHASTIC LOEWNER EVOLUTIONS AND CONFORMAL FIELD THEORIES. 



A. The SLEK/Minimal models connection. 

1. Chordal SLE: definition. 



Throughout this paper, we consider chordal SLE which describes random curves joining two boundary points of a 
connected planar domain. For a comprehensive and detailed introduction to SLE, see e.g. [1, [H, [l^. The definition 
of SLE is most conveniently given in the upper half complex plane H: it describes a fluctuating self-avoiding curve 
7t which emanates from the origin (z = 0) and progresses in a properly chosen time t. If 7f is a simple curve, this 
evolution is defined via the conformal map gt{z) from the domain Ht = IHI/7]o,t], i.e. the upper half plane from which 
the curve is removed, to H. In the more general case of non-simple curves, the function gt{z) produce conformal 
maps from Mt = Kt to M where Kt is the SLE hull at time t. The SLE map gt{z) is uniquely determined by 
fixing three real parameters. A conventional normalization (the so called hydrodynamic normalization) is such that 
gt{z) = z + 2t/ z + ■ ■ ■ near z = co. With this choice the time t corresponds to the upper half-plane capacity CKt of 
the hull Kt, — 2^- The function CK-t is a positive quantity and satisfies the additive law c^^,jg-i^^ ^-j = + ck^, ■ 
The SLE map gtiz) is a solution of the Loewner equation: 

4:9tiz) = —r^ — T .9t=o(2) = ^, (1) 

where is a real valued process, € M, which drives the evolution of the curve. For a system which satisfies the 
Markovian and conformal invariance properties, together with the left-right symmetry, the process is shown [2^ to 
be proportional to a Brownian motion: E[^t] = and E[^t^s] = Kmm{s,t). In the following, we use the symbol E[...] 
for the stochastic average over the Brownian motion. 



2. Martingales in SLE processes. 

In I, i a i 

it was understood how to associate CFT states with the growing curves of the Loewner process. 
To reveal the SLE/CFT connection, the first step is to define stochastic process such that conditioned correlation 
function are martingales for this process. The main idea is to arrange the statistical sum by first summing over all 
the configurations Cj^ presenting a trace 74 of the interface with fixed shape and capacity t/2 and then summing aver 
all the possible shapes of the trace. The first sum defines a conditioned correlation function while the second sum can 
be seen as the stochastic mean of this correlation function. To make it more concrete, consider an observable O of a 
lattice model defined on the domain H. The statistical sum -< O can be written as: 

-< O >n= EH O ^^,] - ^P^'] ^ ^7t: (2) 

7t 
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where -< is the statistical average conditioned to the presence of the trace jt and -P[C^J is the probabiUty of its 
occurrence. The correlator -< O does not depend on the choice of t: the stochastic mean of the correlator -< 
is thus time independent and it is a martingale. 



3. Martingales and conformal correlators. 



The relation ([2]) becomes extremely useful at the critical point where the model is expected to be described in the 
continuum limit by a conformal field theory. In this case, the operator O corresponds generally to a product of, say, 
I primary fields (f>i{zi), 0{{zi}) = Y\\=i 'Pii^i) at positions Zi. The image -fcf) of (j) under a conformal transformation 
f{z) is 



^0, = [^J{z)]^^Mf{z^)), 



(3) 



where A.; is the conformal dimension of . The statistical expectation values can then be expressed in terms of CFT 
correlation fimctions J-{{zi})]si^: 



^{{Zi})Mt = 



< 0({zJ)V^(oo)V^(zt) >Ht 
< -0(oo)V'(zt) >Ht 



(4) 



where < • • • indicates the conformal correlation function computed in the domain = IHI/7t, i.e. the upper 
half plane with the trace removed. The fields ip^Zt) and ip{oo), inserted respectively at the tip zt of 7f and at the 
infinity, are the b.c.c. operators implementing the boundary conditions at which the interface anchor. 

Using the conformal invariance, the correlation functions J-{{zi})ji^ in the domain wall Hi can be expressed as the 
correlation functions in the upper-half plane H: 



< v(oo)V'(6) >H 



(5) 



Note that the Jacobians coming from the conformal transformation on the tp fields cancel between the numerator and 
the denominator in the above expression. 

We are now in the position to understand the SLE/CFT connection. Under the SLE, the trace jt evolves and the 
iterated sequence of infinitesimal conformal mappings gt{z) satisfying Eq.([T]) leads to a Langevin dynamics for the 
conformal correlator J-'{{zi})ji^. Imagine evolving the trace for a time t and then, for an infinitesimal time dt. Using 
Eq.([T]), the variation d{3* (t)^{zi)) =9t+dt (f)^(^Zi) - ^'(f>{z^) is given by: 



= ^ ^ dgt{w)T{w) = 2dt ( - 



A, 



{at{zt)-£,tY gt{zi)-^t 



(6) 



where T(z) is the energy-momentum tensor. The variation for the i/i's are given by the Ito differential: 



dmt)) = d^Mit)dit + ^dli^{^t)dt. (7) 
Using the Eqs.llT])-® in Eq.® and averaging over all the realization of 7f, one obtains the diffusion equation: 



{9t{zi)-£.tY 9t{zi)-£,t 



-di nn{z^}U]■ 



(8) 



As previously said, the correlation function ^{{zi})Mt is a martingale of the SLE process and satisfies therefore the 
following differential equation 



A,: 



d. 



+ H{Z^}U = 0. 



(9) 



4- Operator formalism and minimal models. 

In order to show the consequences of the above relation, it is convenient to introduce the operators L„ which are 
the modes of the energy-momentum tensor T{z): 

LnMz) = J-<{dw w'^'+^TiwUU). (10) 
2t:i ,L, 
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The Ln are the generators of the conformal transformation in the CFT Hilbert space. A primary operator 4>, satisfying 
Eq. ([3]), is annihilated by the positive modes of L„, L„(/) = for n > 0, and the eigenvalue of the zero mode is its 
conformal weight, Locj) = Aip. The L„ satisfy the Virasoro algebra 

[Ln, L,n] = (n - m)Ln+„i + Y^'^^(" - l)'5„+m- (11) 

All local fields of the theory can be obtained by applying the L_„, n > 1, to a primary fields <j). Each set of states of 
the form L^j^ • • • Lk„^ 4> forms a conformal family [0] and corresponds to the representation of the highest weight of the 
Virasoro algebra, also called Verma module, the primary field corresponding to the highest vector. The descendant 
state L-ki ■ • • ^-fc,„0: X] ~ s&^A to be the n—ih level of the module of and its conformal weight is A + n. 
For c < 1, the unitary representations of the Virasoro algebra has central charge [2H : 

c=l-—^— m-2,3,--- (12) 
mym + Ij 

and correspond to the simplest family of conformal theories, called minimal models Mm . The minimal model Mm is 
characterized by a finite number of Verma modules [<j)r,s\ with conformal weight: 

[(m + l)r — msl^ — 1 , , 

^r,s=^^ , ' , r = l,2,...TO-l, s = l,2,...r. 13 

4:m[m + 1) 



5. Null vectors in CFT and SLE. 



In CFT, the differential equations satisfied by the correlation function can be derived from the structure of the 
representation modules and, in particular, from the relations between the descendant states. One can show that 
^{{zi})wt is a solution of ([9]) if the b.c.c. operator ip obeys the condition: 

(i_2 - ^L\)i; - 0, (14) 

The equation Eq. (fT4|) means that the descendant states L'^i4' and L_2'0 are linearly dependent. The operator ip 
transforms then as a degenerate representation of the Virasoro algebra which has a null state at level 2. The operator 
Tp can thus be identified with the operators (^1^2 or (f)2^i of the minimal models Mm which are shown to satisfy: 

(L-. - ^(^^^L^_.)0r. = 0, (15) 

(L-. - ^(^^^i^^)^.! - 0. (16) 

Comparing Eq.ljTS]) and Eq. ^TB]) to Eg. p^ and using Eg. lT^ and Eq. fT^l) . one can finally make explicit that the 
SLE/CFT connection: an SLE process with parameter k describes interfaces in CFT with central charge: 

(6-^)(3«-8) 

^ 2^ ' ^'^^ 

and the boundary conformal operator ip has scaling dimension: 

A. = (18) 
2k 

This derivation shows that martingales for SLE processes are closely related with the existence of null vectors in the 
the appropriate Verma module. 



III. SLES FOR CFTS WITH EXTRA LIE GROUP SYMMETRIES. 



As mentioned above, the operator content of the minimal models Mm, which have central charge c < 1, corresponds 
to the space of representations of the Virasoro algebra (jlip . In this case, the conformal symmetry determines 
completely the spectrum of the theory and the differential equations satisfied by the correlation functions. 

The applications of the SLEs defined in the previous section are limited to CFTs with central charge c < 1. This 
can be directly understood from the fact that, for a CFT with extended symmetry, not all the local fields can be 
obtained by the application of the L„. In particular, one expects at each level additional states and the linear relations 
of the type P^ - P^ will include additional terms. 
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A. SLEs and CFTs with Lie group symmetries. 



The SU (2) WZW models are among the most huportant CFTs with additional symmetries. These CFTs posses an 
internal SU (2) continuous symmetry which, together with the conformal one, is realized by the set of chiral currents 
J°'{z), a = +,—,0, with scaling dimension Aj = 1. The current algebra, derived from the current-current operator 
product expansion (OPE), takes the form: 

[Jn. JL] = fc^'Jn+m + Nn5„+m, (19) 

where are the current modes, — 1/ {2Tii) §_^dww^ J°'{w), the are the structure constants of the su{2) Lie 
algebra and N is the level of the algebra. The representations of the Kac-Moody algebra , leads to a family of CFTs, 
denoted SU{2)n, with central charge cn = 3N/{N + 2). The local fields (l)j{z) of the theory transform under the 
action of the chiral currents in the representation space of dimension j{j + 1) of SU{2). Analogously to the Virasoro 
case, a WZW primary field is defined by: 

j;j0,(z) = O forn>0; JoV,(z) = <>,(z), (20) 

where the i° are the SU{2) generator matrices in the j representation. The conformal weight Aj of a primary field 

(f>j{z) is 



Lo0,(^) = A,^,{z) = ^^±Il0^.(z). (21) 



Under an infinitesimal conformal dg{z) transformation together with an infinitesimal SU{2) gauge transformation 
d9°'{z), the primary fields (pjiz) changes as: 

d{(l)j){z) = — (f dw dg{w)T{w) (f>,{z) + — I dw dB°-{w)r{w) <pAz) (22) 
2m 2m 

In [l3|, the SLE approach has been extended to these theories. It proposes to describe this model by a composition 
of two independent Brownian motions, one in the physical space, defined as in Eq.lll), and one in the internal SU (2) 
spin degrees of freedom: 

dOtiz) ^ /f , n0'l0l] = r5'^^'unn{s,t), (23) 
9t{z) -£,t 

where the variance t is an independent parameter of the SLE evolution. The idea is then to attach to trace 74 a spin 
1/2 degrees of freedom. The interface then undergoes both a standard SLE evolution in the physical space and a 
stochastic SU{2) rotation. 

The derivation of the SLE/CFT connection is strictly analogous to the one shown before. The difference here is 
that the b.c.c. operator V'(Ct) Eq.Q carries a spin 1/2, tp = V'i/2(^t) which transforms under the representation 
4>i/2- Using Eq.([T|) and Eq. (I23p in Eq. ip^ . the Ito formula for ■01/2 (Ct); written in terms of the modes L„ and J^, is: 

d{i^i/2m = ^6^-1^1/2(6) + d^V^ 17/^1/2 (6) + dt (^^L\ + ^ J^i^^i) V^i/2(6)- (24) 

The operator O is now a product of WZW primaries and remain fixed while the trace evolves. This means one has 
only to consider their variation ([6]) under the conformal transformation gt which maps the domain Ht in the domain 
Hj. The martingale condition (fT4|) takes the form: 

r 



i-2--L^ 



'-J-ir-i) ^112 = 0. (25) 



The correlation functions in WTjW theories satisfy first order differential equations, called Knizhnik-Zamolodchikov 
equations, which are derived from a relation between descendant states at the first level: 

L_i + .^^J-^e^ ^ 0. (26) 

Applying the operator L+i to Eq. ([25|) and using the commutation relation ([19)) together with [Ln J° ] = — mJ^_,_„, 
one fixes a first relation between k, t and N, t/{6 — k(2Ai/2 + 1)) = 2/{N + 2). It was shown [17j that a second 
condition, k + t = 4, can be obtained by demanding the one-point function of the current to exhibit a simple pole at 
the tip of the trace zt- The final result is: 
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IV. SLES IN Z{N) PARAFERMIONIC THEORIES. 

We consider the CFTs with extended Zjv symmetries, the so calles Z{N) parafcrmionic theories. In particular, 
we focus on the Z{N) parafcrmionic theories with central charge c = 2{N — 1)/{N + 2) (c > 1 for > 4). These 
theories describe the continuum limit of the A'^-states spin models interacting via a Zjv invariant nearest-neighbor 
coupling, at their self-dual critical points (see section |V]). For N = 2 and iV = 3 one finds the well known Ising and 
three-state Potts model. The corresponding CFTs has central charge c = 1/2 and c = 4/5 respectively and coincide 
with the minimal models M3 and M5 : the operator content of these theories can also be determined by studying the 
representations of the Virasoro algebra. In this sense, the Z2 and symmetry of these models are trivially realized. 
The CFTs describing the critical point of the Ising and the three-state Potts are then expected to be described by 
the SLE in Eq.([T]). The b.c.c. conformal operators implementing the SLE interface boundary conditions have been 
identified. We will show that these results admit a natural extension to general N. 



A. Parafermionic current algebra. 

We briefly review the Z{N) parafermionic theories with central charge c — 2{N — i)/{N + 2). These theories 
were introduced and constructed in (lol . [2^ . Here we shall enounce the main results in a slightly different manner. 
The general arguments and the notations used here are strictly analogous to the ones discussed in a series of papers 
[H, HJ, m, nil where a second series of parafermionic theories was studied. 

Extra Zn group symmetries in two-dimensional conformal field theories are generally realized by a set of holomorphic 
currents '^'^(z) (fc = 1, • • • , — 1), satisfying the following operator product expansion: 

\k.k' 

(z-z')^+ 
X {*'^+'='(z') + 0(z-z')}, k + 

]_ {. , ,.o2A* 

{z-z'f 



^\z)^-\z') = — ^ ^ + (z - z' f^T{z') + . . . ^ (29) 



where A* is the conformal dimension of the parafermions "^k{z) and A^^,;,, are the structure constants of the algebra. 
We shall be interested in parafermionic theories in which the dimensions A* of the parafermions {5''^} take the 
minimal possible values admitted by the associativity constraint: 

A* = A*_, = M^^ fc = 0,l...iV-l. (30) 

Note that the above formula is not symmetric under the exchange fc — fc : the field '^^^ in (p9)l is assumed to have 
dimension A* in the sense that the indices fc referring to the Z^f charge are always defined modulo iV. Thus, 

^ ^ (*'r , A*_,^A*,. (31) 

The structure constants A^i^^., and the central charge c (of the Virasoro algebra) are given by the expressions: 



,.Kk' .2 _ (fc + fcO!(iV-fc)!(iV-fcO! 

^ ^ " fc!fc'!(iV-fc-fc')W! ^ ^ 

2(^-1) , . 



B. Representation space: the Zjv sector. 

The CFTs we are considering describes the self-dual (Kramers- Wannier invariance) critical point of Zjv-invariant 
lattice model 10|]. In addition to the Z^r symmetry, the theory possesses the dual Zj^ invariance, which is enanched 
by the conservation of the antiholomorphic fields ^{z\ The Hilbert space of a Z^^ x Z-^ invariant theory, splits 
into subspaces characterized by the Zj^ x Zjv charges {p, g}. A field (^{p.g} in one of these subspaces transform as 
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(j)^pqy{z) exp[2iTr{pin+qn)/N]4>p^q{z) under global rotations of angles 27rm/iV and 27rn/iV". The lattice spin operator 
Cfe (see|Vl) and its dual /i^ are described in the continuum limit by fields of charge {k, 0} and {0, k} respectively. The 
(anti-)holomorphic currents ^''^(z) ( ^'(z)), which appear in the OPE of atlJ-k (cfc/i^), have charges {k, k} {{k, —k}). 
Recently, the lattice holomorphic realization of the parafermionic currents has been discussed in [3]|. 

It is in general convenient to express the charge fc, k' as [q* ,p*] = {k -\- k' ,k — fc'}, where now q and p are defined 
mod 2N and q* + p* even. In this notation, the currents 5''^ and ^'(z) have charge [2fc,0] and [0, 2fc] respectively. 

In the following we concentrate on the action of the holomorphic field ^''^(z), all the results being valid also for ^'('2). 
Thus, without losing any generality, one can study the structure of the Zn representations which we denote with 
<I>^' . Moreover, we find convenient to use a different convention q for the Zn charges, defined by q* — 2q mod2A^. 
With this choice, the product is a field with charge k ± q. 

For N odd, we consider then the representation fields 

g = 0,±I,--- ,±(iV-l)/2 (TV odd). (34) 

For TV even, it turns out that the modules of the representation corresponding to ^^"^ with [Af/4J < q < [N/2\ are 
identical to those of < g < [A''/4J. In order to recover the right number of Zn representations one has to consider 
half-integer value of q: 

<i>±«(z,z) <7 = 0,±i,±l,--- ,±LiV/4j (TV even). (35) 
Naturally, the physical meaning of the the half-integer charge is recovered in the usual notation q* . 



1. Currents modes in Zn sector. 

The currents {^''^} can be decomposed into mode operators v4*:.^„, whose action is to change the Zn charge of the 
representation fields 

SI = niod 1. (37) 

The value is the first level in the module of corresponding to the Z^ charge k: it determines thus the 
level structure of the modules induced by the Z^ representation fields. The values of $1 can be easily obtained 
by considering the module of the identity whose descendant states are the chiral currents '^''{z). The levels of 
these operators correspond to their conformal dimensions A*^,. Taking into account that, owing to the Abelian 
monodromy of the fields ^''^(z) in the Z^ sector, the level spacing is equal to 1, one can readily obtain from Eg. ([50]) 
(5° = — fc^/TVmod 1. The level structure of a generic field $='=*(z), Ea. ([57)) . is then extracted from the module of the 
identity by inspecting its corresponding submodule. Note also that = for = —2q. This results in the presence 

of zero modes A^"^: 



AJ29$±9(0) = /ig$«(o). (38) 

which associate at each field with conformal dimension A*, the field <f>^'' with opposite charge and with the 
same conformal dimension. A* = A^^. The eigenvalues hg defined in (j38[) characterize the representations of the 
parafermionic algebra together with the conformal dimension the fields As usual, primary fields are defined by 
^'Lxi I = for n > 0. Each representation module is then characterized by the two primary fields <1>^* which 

can be obtained one from the other as in (|38p . An example of the module structure associated to the field <I>^^ and 
of the action of the currents is given in Fig[5]for N — 5. Using the formula ((37)) . we have 6^^ — 1/5 and = 4/5. 
Using the expansion (|36p . the action of the modes in each sector can be given in terms of a contour integral: 

Atsi^ ^^f"^' (z)^*-*-'+.+"-ivi;'=(^)$«(0). (39) 

In the Zff sector, one can simply consider the action of the currents {^'^^} because, via the Eq. (I^S)) . they completely 
determine the full algebra of the other currents {'^'^''}, fc = 2, • • • , [A^/2J . 
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The commutation relations of the mode operators can be deduced from Eg. ((39)) by using standard techniques in 
the complex plane. These relations are given in Appendix IB II 

As shown in the Appendix IB H the representation space contains [iV/2j+l primary operators with q = 
0, 1, • • • , f for iV odd and with q ^ 0, ±i, ±1, • • • ,±IN/A\ for N even. The conformal dimension of these oper- 



ators turn out to be: 



^^= N^N + 2) 1 = 0,l,---.^iNodd), 5^0,1 (iVeven). (40) 



As said above, to each primary $^ corresponds another primary field <i>^'? of the same dimension. 

In the case N = 2 and = 3 (i.e. Ising and three state Potts model), one has only two primaries of the 
parafermionic algebra. For N = 2 one finds the identity operator with Aq = and the operator ^^^^"^^ with 
dimension A±i/2 = 1/16, corresponding to the Ising model spin operator. Analogously, for A^ = 3 one has, together 
with the identity operator, the operator with A±i = 1/15. As previously said, the case N — 2 and A^ = 3 
are special because they can be identified with the minimal model A/a and M^. The operators ^^^/^ and are 
identified respectively to the operators (/)i_2 and (j)2,3 of the corresponding M3 and M5 Kac tables, see Ea. (|T3l) . As 
we will show more in detail below, the representation module corresponding to contains one descendant state, 
£(-'v=3) _ ^~i^^<j)+i = A^_^^^^~^ singlet under Z3 transformation (i.e. q = 0). This singlet (energy) operator has 

dimension Ag = 1/15 + 1/3 = 2/5, it is a primary of the Virasoro algebra and can be identified with the operator 

£(Af=3) = 02 1 of the M5 Kac table. 



C. Representation space: the disorder sector. 



A key observation is that the theory we are considering is actually invariant under the dihedral group _Djv which 
includes Z]y as a subgroup. This can directly be seen from the symmetry of Ea. (P5)l - ([5^ under the conjugation of the 
Zjv charge, q ^ N — q. At the level of the lattice model (see|V|, this comes from the invariance of the Hamiltonian 
under the transformation <Jk{x) — > (T^(a;). 

The space of representation thus includes the A^— plet of Z2 operators which we denote as: 

{Ra{z,z), a = l,...,N}. (41) 

In the following, we refer to these operators as disorder operators. 

The theory of disorder operators has been fully developed in I22|. A detailed discussion about the general properties 
(products, analytic continuations) of these operators is given in [23j . 

The most important feature of the disorder operators is their non-abelian monodromy with respect to the chiral 
fields {4'='='^}. This amounts to the decomposition of the local products \E''^(z)i?a(0) into half-integer powers of z: 

■^\z)R,{0) = Y,—^A\RM. k=l,2,...,[Nl2\. (42) 
[ z ] ^ 2 

The expansion of the product ^~^{z)Ra{Q) (with fc = 1,2, . . . , [A^/2J) can be obtained by an analytic continuation 
of z around on both sides of Eq. (|42p . The result is: 

^-\z)RM ^ , J^r. 4 Ra^m. A; = 1,2,...,LA^/2J. (43) 

^(z)^fe + 2 

In accordance with these expansions, the mode operators can be defined by the contour integrals 

^ii?a(0)--^ / Az{zt"^^--^^\z)RM. (44) 

where the integrations are defined by letting z turn twice around the operator i?a(0) at the origin, exactly as described 
in (2^ . [2^ . The commutation relations between the current modes in this sector can be computed and the dimension 
of the primary disorder operators is (see Appendix IB 4^ : 

Af = ^^^^#^^ . = 0,l,...,LAr/2j. (45) 
16(Ar + 2) ' . . L / J \ ) 



9 



Note that iov N = 2 one has the operators and R^^ with dimension = 1/16 and Af = 0. This is of course 
expected as the dihedral group D2 = Z2, i.e. the cychc Zjyf elements and the reflection Z2 elements of the group 
coincide for N = 2. In particular one has i?° = $='=1/^ while coincides with the identity. For = 3, the disorder 
sector contains the operator {Aq =1/8) and R^ {Af = 1/40). In terms of primaries of the Virasoro algebra they 
coincide respectively with the 0i^2 and (j)2,2 operators of the M5 Kac table, see again Ea. p3)) . 



D. Stochastic motion in the internal space. 



We discuss now a possible SLE approach to describe parafermionic theories. In the previous section we have seen 
that the action p9|l and (|44|) of the chiral currents modes on the representation fields amounts to a twist either of the 
Zjv charge q for fields ^^(z) belonging to the Zjv sector, or of the Z2 index a for fields Ra{z) in the disorder sector. 
The b.c.c. operator is expected to transform as an operator in the Z^ or in the Z2 sector. By analogy to the case 
of the SU{2)n WZW [17|, we consider an independent stochastic motion in the additional degrees of freedom. We 
provide a parafermionic field the following evolution: 

X{z) X(z) + / dgtiw)T{w)X{z) + / dwidOf (w)-^^ {w)X{z), (46) 
2m Jq zttz Jq 

where X{z) = or X{z) = Ra{z) and the contour Cz is a general closed contour around z (remember that in the R 
sector the Riemann surface is two- fold). This evolution contains, in addition to the variation under an infinitesimal 
conformal transformation dgtiw), a component with shifted Z^ charge q q + k ot Z2 charge, a — > a + /c. The (j46[) 
can be seen as a motion in the parafermionic representation module. In the sequel, we shall concentrate on considering 
only the action of the most fundamental 'I'*^ fields. 

The definition of the internal stochastic motion is given in such a way that the corresponding Ito derivative of 
ipiz) takes the form ([7]) or (|24p . More precisely, we want to derive a martingale condition as a linear relation between 
descendants at the level two of some representation module. This allows for a direct comparison with the CFT results. 

In the case the boundary operator ip{z) coincides with a primary in the Zj^ sector, we define: 

(^^)* (^) = . . . !1^A^+6-' +1 ' ^^^'^s'] ^rS^^-'^' min(.,t). (47) 

The exponent in the denominator has been chosen according to Eq. (j36p which encodes the OPE between the currents 
and the fields Note that this evolution is very similar to in the su{2) Lie algebra. This is somehow 
expected since the Z{N) parafermionic theory corresponds to the coset SU{2)ff/U{l) (see appendix E]) . 
In the case when ip{z) transforms as a disorder Ra field, we define, according to Eq. fi^ : 

(^^)*(^) = rTT^^7V^^' netO':']^r5^--'^'nMs,t). (48) 

We point out that the definition of the stochastic motion given above is by now purely algebraic and based on a 
formal analogy with the case SU{2). However, the CFT predictions and the identification of the possible interface on 
the lattice model motivate the interest of such construction. 



E. Martingale conditions and parafermionic results. 

As we will discuss in some detail later, the Ising model {N = 2 and c = 1/2) shows in its FK representation (see lyj) 
a domain wall described by the SIjEiiq/^. The b.c.c operator creating such interface is identified with the operator 0i^2 
of the M3 Kac table, what is consistent with the Eq. ([TT]) and Eq. ([TB]) . In the parafermionic representations, we have 
seen that this operator corresponds to the disorder operator or, equivantely, to the $±1/2 operator (remember 
that iV = 2 is special as D2=Z2). In the three state Potts model {N = 3), the operator (j)i^2 of the M5 Kac table is 
again identified again with the disorder and it is expected to create an SLE24/5, in agreement again with Eq. (jl7p 
and Ea. p8)) . The operator 02,1, on the other hand, can be associated to an interface described by the SLE10/3 and 
it corresponds to the e'^^^'^^ — A^_^^^^~^ — A~^^^^^ operator mentioned before. 

On the basis of these correspondences, it is then natural to identify for general N the b.c.c. operator ip{z) with the 
disorder primary operator with dimension Aq = {N'^ + N — 2)/(16(A^ + 2)) and with the Zjv singlet (i.e. charge 
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(7 = 0) operator e*^^) = A_^^j^^^ ~ A^_^ij^^ ^ with dimension = IjiN + 2) . Note that this singlet operator is 

not a primary of the parafermionic algebra but it appears as descendant in the module of <i>*^ (see appendix [B|) . 
Using Eg. ([50]) . Eq. ([57|) and Eg. ([^5]) . the stochastic motions ([TT]) and amounts respectively to the following 

dynamics for "^(z) in the two cases when "^(z) = ^e= find "tpiz) = ijjR = R^'a^'- 



(49) 



and ip{z) = -tjjj^ = i?!"-*: 

rf(^i?(6)) = d^tL-M^t) + {d9)^A\,ijR{^t) + {d9)-^AzliJB.{^t) 
+dt {'^Ll, + ^ [AI.AZI + AzlAl,]) ^«(6). 

The martingale condition, to be compared with Eg. (fT4|) and Eq. ljSS]) . now reads: 



( 



"^-l/jV-l"^l/Af-l 



(50) 

(51) 
(52) 

We can now enounce the main results of this paper. We have studied the representation modules of the e and i?!"' 
operator. For iV > 4 one finds that there are at least two independent states at the second level of these operators. 
We have found the following relations at the second level of these modules: 



L-2 - -^Ll, - ^ [Al^A-_\+Az\Al^]) M^t) 



N+1 
N + 2' 



A 



i-2 - 



iV + 2 



l/N-l^'-l/N 
2-4/JV^3 







A{N + 2){N - 



-^[Al,Azl+AzlAl,]^^Ri^,) = 0, 



(53) 
(54) 



The derivation of these relations is shown in appendix |B] We have verified that, for N = 2, the relation ([M]) reduces 
to p3|) and for = 3, the Eq. ljSS]) and Eg. ([54]) are respectively equivalent to Eg. lfTH)) and Eq. lfTS)) . In these cases, 
which as already said coincide with minimal models M3 and , all the descendants of the theory can be obtained by 
applying negative Virasoro modes. In our case, this amounts to the fact that the states at the second level constructed 
from the application of parafermionic modes can be written as a linear combination of and L_2. 

Comparing Ea. ([5T|) and Ea. ([52| to Ea. ([53| and Eq. ([53)) . we obtain, for > 4 and for the two SLEs evolutions: 



N+1 

A— — - Tl = 4 



= 4 



iV- 
.N - 



{N 



N +1 



T2 = 2-4/^ 



[N + 2){N + l) 



(55) 
(56) 



V. Z{N) SELF-DUAL CRITICAL SPIN MODELS. 

Consider a square lattice with the spin variables (jj at each sites j taking N possible values: 



CTj — exp 



i27r ■ 



n(j) = 0,l,--- ,iV-l. 



(57) 



The mo st g eneral ZAr-invariant spin model with nearest-neighbor interactions is defined by the reduced Hamiltonian 



VN/2\ 



/ 2'Kmn\ 


- 1 


hi N ) 







(58) 
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and the associated partition function reads: 



cxp 



-(3 H[n{i) - n{j)] 



(59) 



In fact, the Potts model is recovered in the case — J for all m. The model presents, in this case, a permutational 
^jv symmetry. Another meaningful model is the clock model which is defined by the partition function (|59p with 
Jm ~ J5m,i- Defining the Boltzmann weights: 



exp [~(iH{n)] , 



0,1, ••• ,iV-l, 



(60) 



the general Z{N) spin system is described by [-/V/2J independent parameters x. The Kramers and Wannier (order- 
disorder) duality has proven to be a powerful tool for examining the behavior of these models. The phase diagram of 
the model ([55]) for = 5 states of spin is shown in Fig. ([2]). It has been show that the Z{b) model is self-dual on the 
line xi+ X2 — 1/2(V5 — 1). Here we are interested in the existence of two critical points (x*, 0:2) and (0:2, x*), where 
x\ = sin(7r/20)/ sin(37r/20) « 0.34 and X2 = sin(7r/4)/ sin(77r/20)2:i « 0.27, where the model is completely integrable. 
The Z(b) parafermionic theory defined in the previous sections describes the continuum limit of the lattice model at 
these points. In general, in the self-dual subspace of which contains also the Potts model the clock model, 

the Z{N) spin model is completely integrable at the points [ll| : 



Xq = 1 



fc=0 



n/ Trfc TT 



sm 



7r(fc-M) 

N 



47V 



(61) 



These critical points are described by the Z{N) parafermionic field theories. 



A. Boundary conditions and SLE interfaces. 

Consider the three-states Potts model. From the Eq. pT]) and the Eq.((T8]), as previously mentioned, the b.c.c. 
operator T/ie^""^-* is expected to create an SLExo/3 trace. In this paragraph we suggest a definition of the interfaces 
within the lattice models whose scaling limit is described in these SLE processes. In the following we indicate the 
possible values of the spins with the letters A,B---. The b.c.c. operator has been shown to generate the 

boundary condition where the spins are fixed to (say) the value A on the left side of the origin while can take the 
value B or C with equal probability on the right side [l3i |. With this boundary condition, to which we refer with the 
short-hand notation A\B + there exists a single domain wall between the spin A and the spins B and C, as shown 
in FigUJ The statistical properties of this interface has recently been analyzed by numerical means in |32l | and are 
shown to be well described in the continuum limit by an SLE10/3. 

The generalization to the iV-states spin model is straightforward. Indeed, one can show (see appendix [C| that 

the b.c.c. operator V4^'' produces the boundary conditions A\B + C + £> + ••• , where the spins on the left of the 
origin take the value A while the ones on the right take with equal probability the — 1 values of spin different from 
A. As in the case of A^ = 3, there exists a single boundary domain between the A spin and the B,C, - ■ ■ ... Note 
that the boundary conditions do not satisfy the reflection symmetry and the measure on the curve should reflect this 
asymmetry. However, in the continuum limit the reflection symmetry is expected to be restored as explicitely shown 
in [3^. This can be argued from the fact that, in this limit, the statistical expectation values can be expressed in 
terms of CFT correlation functions which contain the b.c.c. operators, see Eq.(|4]). These correlation functions satisfy 
indeed the reflection symmetry. 

It is interesting to note that one can assign a "color" to the interface depending on the values of the spins {B,C,- ■ ■) 
at the right of the domain wall. The point is that, for A^ < 4, the Boltzmann weights a;* , n = 1, 2, • • • , A^ — 1, are 
equal: the energy of one bond does not depend on the difference between the values of the nearest-neighbor spins. In 
this sense, the Z^ symmetry is trivially realized and this is essentially the reason why, for A^ = 2, 3, the correspondent 
CFT can be constructed from the conformal invariance alone. However this is not true for A^ > 4 where the Boltzmann 
weights X* , n = 1, 2, • • • [A'^/2J are different. In the continuum limit, this amounts to a more structured chiral algebra 
with a central charge c > 1. From the point of view of the lattice model, one can observe that for A^ = 2, 3 the energy 
cost to have a spin i? or C at one side of the interface is the same. Thus, as A^ > 4, one expects the internal spin degree 
of freedom to play an important role in the spatial evolution of the interface. This is somehow in agreement with the 
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) o o 



• o • • • • 



b) 



FIG. 1: a)Three state Potts model: interface (green curve) between the A spins (black circles) and the B (full blue circles) and 
C (empty blue circles). The boundary conditions A\B + C as in the text, b) The interface defined in the Z{5) spin model with 
boundary conditions A\B + C + D + E. 



need to introduce an additional stochastic motion in the Zn group. However a clear the geometrical interpretation 
of the stochastic evolution in the representation modules (see section H VP I) is missing. 

From the results (|53|) and ((55|l , we hypothesize that the geometric properties of this interface for iV > 4 are described 
by a SLE with k = A{N + 1)/{N + 2). In particular, since k < 4, we expect the domain under consideration to be a 
simple curve for each N. 

We discuss now the identification of the interface on the lattice associated to the b.c.c operator tpji, see Eg. ([54)1 . 
In the Ising and three states Potts model, this operator is associated to an and SLE24/5 respectively. The 

correspondent interfaces on the lattice can be identified with the help of the Fortuin-Kasteleyn (FK) representations. 

The partition function for the N = 2 and TV = 3 critical models take the form: 



^exp 



-13 J2 H(n{€) - n{j)) 



<ij> 



XI n [^1 + (1 - a;*)(5„(i),„(j-)] , 
{<t} <ii> 



(62) 



where = \/2 — 1 or = x*_^ — (VS 
respectively. By expanding the product 
cluster representation: 



l)/2 is the critical point of the Ising and of the three states Potts model 
and summing over all the spin configurations, one obtains the FK random 



Z 



c 



^\M-b/ 



(63) 



where the sum is over the subgraphs C oi G, C C the graph Q = {V,E) being composed by the V sites and the 
E edges of the square lattice. Each graph C is specified by the b bonds on the edges and the c clusters (=connected 
component which include the single sites) in a given bond configuration. 

Consider the following boundary conditions: all the edges on the negative axis carry bonds while all the edges on 
the positive axis carry no bonds. For the spin variables this means that all the spins located at the left of the origin 
take the same value while the spins on right side are unconstrained, i.e. they are free (F) to take all the N possible 
values. 
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FIG. 2: Phase diagram of th Ziji) spin model 

In such an arrangement, each graph C presents a cluster growing from the negative axis into the upper half-plane. 
The boundary of this cluster, starting from the origin, is conjectured (see for istance [ll], chapter 5) to be statistically 
equivalent in the continuum limit to the SLEig/3 for N = 2 (Ising model) and to the SLE24/5 for = 3 (the three 

states Potts model). This is in agreement with the properties of the b.c.c. ij}^^^ operators, N — 2,3,4, which have 
been shown [s^ to be associated to the A\F boundary conditions. 

The generalization to the case iV > 4 is more cumbersome. To our knowledge, the modular properties of the 
i?— sector characters, and thus the boundary conditions associated to these operators, have not been yet studied for 
these cases. This motivates a more detailed analysis of the boundary conformal properties of this sector. However, 
we point out that at the level of the parafermionic algebra representations the main features of the ij^R fields, such as 
the structure and the degeneracies of the representation module or the current zero modes eigenvalues (see Appendix 
IB 4p . directly generalize to each N > 2. The formula (|45p is an example. It is thus quite natural to suppose that the 

operators ip^j^"^ produce the A\F boundary conditions for each N. Another observation is that it is quite direct to find 
an alternative formulation of (j59p in terms of colored clusters on a graph. For sake of clarity, we focus our attention 
on the case N = 5. 

The partition function for = 5 at the critical point a;* , see Eq. (pT|) . can be rewritten as: 

^(JV=5) ^Y^li [a^t + (1 - ^l)Sni.)M,) + i^*2 ~ ^l)Sni^)Mn±2] ■ (64) 

{cr} <y> 

If xl — X2 one recovers the five states Potts model in the usual FK representation. 

Expanding the factors in ([M)) . one has to take into account two type of bonds: the bonds (type 1) connecting spins 
which take equal values and the bonds (type 2) connecting spins which differ by an angle 47r/5. Then, summing over 
the spin configurations, the partition function can be cast in terms of a sum over "colored" graphs C: 

ZiN=5) ^ J2{j:*i)^-bi-b2(^i _ xlf^ix* ~ (65) 
G' 

Each graph C is specified by the number c of connected component and by the numbers bi and 62 of the bonds of 
type 1 and type 2. Naturally, the sum over the differences of the spins along a loop must be zero (modulo 27r). For 
instance, a loop formed by three bonds of type 2 and one bond of type 1 is not permitted. The prime in C indicates 
that the sum in (j65p is taken over the allowed graphs. A typical configuration is shown in Fig. ([3]). 

By analogy to the Ising and three states Potts model, we consider the boundary conditions in which the spins have 
a fixed value on the negative axis while they are free on the positive axis. With such boundary conditions, there is a 
cluster of bonds of type 1 growing from the negative axis, see Fig. ^ . 

From (j56p . we conjecture the geometric properties of the boundary of this cluster are described in the continuum 
limit by the SLE14/3. 
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FIG. 3: Cluster expansion of the Z{5) spin model. The spins at the vertices of a black bond have the same value while the 
ones at the vertices of a red bond differ by an angle 47r/5 




a) 



FIG. 4: a) A typical cluster configuration of the Z{5) spin model with type 1 (black) and type 2 (red) bonds. The boundary 
conditions on the real axes are as in the text. The green curve is the boundary of the cluster growing from the negative axes, 
b) A piece of the interface which crosses red or empty bond. 



The above discussion can be directly generalized to each value of A''. In particular, the Z{N) spin model admits a 
random cluster representations with [{N — 1)/2J types of bonds. Imposing the conditions on the real axis discussed 
above, we conjecture that the measure of the boundary of the cluster growing from the negative real axis is described 
by an SLE with k = 4(iV + 2)/{N + 1). 

Finally we remark that, during its evolution, the interface crosses edges without bonds or with bonds of the type 
2, 3, • • • [{N — 1)/2J which have different energetic costs. Analogously to the case of the interface generated by the 
A\B + C + ■ ■ ■ boundary conditions, one can color the interface depending on the type of bond it crosses. Again, 
the interface is naturally provided of an Zn additional internal degree of freedom, in agreement with the approach 
proposed in Section (|IVD[) . 



VI. CONCLUSIONS. 



In this paper we have considered the possibility to extend an SLE approach to CFTs with additional Zn discrete 
symmetries. These theories describe the continuum limit of the self-dual critical Z{N) spin models. The case N = 2 
and iV = 3 correspond to the Ising and three states Potts model where the identification of the interface on the 
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lattice and of the correspondent b.c.c. operators are known. Using these results, we have identified the possible b.c.c. 
operators for general N. We show that these operators satisfy a two level null vector condition. For > 4, an 
additional term enter these relations. This term is obtained from the action of the parafermionic currents and, for 
> 4 (c > 1), it cannot be derived by the Virasoro modes action. These results led us to propose an additional 
stochastic motion in the different charge sectors and in the Z2 sectors which describe the representation modules 
of the parafermionic algebra. We have finally discussed the possible interfaces in the lattice to which the CFT/SLE 
results should apply. The fact that SLE candidates have been identified in lattice spin models described by non- 
conformal minimal theories opens the possibility to study the SLE/non-niinimal CFT connection by a numerical 
approach. 

It will be interesting to verify our predictions by numerical studies. Further we believe these studies can also better 
clarify the geometrical interpretation of the additional stochastic motion in the internal space. 
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APPENDIX A: RELATION BETWEEN THE SU{2)n WZW MODEL AND THE Z{N) PARAFERMIONIC 

THEORY. 



In this appendix we briefly show how the stochastic motion (pS)) for the SU{2)n WZW theories can be somehow 
compared to the one ((47|) by taking into account that the Z{N) parafermionic theory corresponds to the coset 
SU{2)i\f /U{1). The S'C/(2)jv WZW theory can then be seen as a combined system of the Zjv parafermionic theory 
and one free massless boson field (t){z). The chiral currents generating the Kac- Moody algebra (|19l) can be written in 
terms of the parafermionic currents and the free field as: 



J± = 



A^*^^(z)exp 



±i 



1 



J" 



'N 



(Al) 

yC/(l) 



The stress-energy tensor T-^^^^)" and T^^^^ of the two theories are symply related by T^u(2)m ^ j^z(n) 
where T^^^) = ~l/'i{d(l)Y is the stress -energy tensor of the free bosonic theory. 

All the primary fields (pj of the WZW theory can be obtained by a product of parafermionic fields in the Zm sector 
and of vertex operators of the bosonic theory: 



(l)j{z) ^ (f>^ exp 



J 



N 



(^) 



(A2) 



Using the fact that the conformal dimension 



can easily verify that Aj 



I \vert. 



see Eq. 



of the vertex operator exp[±iQ;(/)(2;)/VA'] is /Sy^-^*-- — a^/N, one 
2T|) and Eq (gD]). Note that the field (jjjiz) in Eq.jASJ satisfies 



jQ0j = jcj)™. In the following we denote (j)J^{z), with m = —j, + ,j, the field in the SU{2) representation of 
spin j such that Jqc/)"^ — m 



vertex operators as: 



The field 



can be expressed in terms of the parafermionic primaries and of 



= a>J")exp 



(^) 



(A3) 



where 61 has been defined in Eg. ((37)) . For instance, the SU{2) operator ipl corresponds to e(^\ see section |lVll 
while the operator (p^^ to exp[±i(j){z)] 

Using the above relations, one can compare the definitions (j23p and (I47p . and in particular the exponent in the 
denominator of the two formulas. This exponent actually determines the current mode acting on the b.c.c. field under 
his evolution. In the SU{2)n theory this is the J^i mode, as it can be seen in Eq.([M|). Let us consider for istance 
the state JZ 



j: 



10} = l/{2ni) (f dziz-wy^ J-{z)(Pl{w) 



J 1L 



l/{2m) j> dz{z - w)-^ (*"i(z)$i(w)) ^exp 



(X l/{2m) f dz{z - u;)-i-2/iv (^-i(^)$i(u;)) ( 1 + _ w)d^{w) + 



—I 

7a' 

i 



(^) 



exp 



(A4) 
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where we have used the relations 
exp " 



2j) and the divergences coming from the OPE in the bosonic sector: 



N 



exp 



{z-z') 



exp 



.q + q' 



0(z') 



(A5) 



The above expression gives precisely the state Aj^jj^ .^^^^ multiplied by a bosonic piece. In the same way one 



can show that the term J_i4>l and the term {J'^^J_^ + J -^ jj^-^)^" correspond respectively to A^^^^^e'^'^^ and to 
(^ii/AT-i^i/JV-i ^-i/AT-i^i/A'-i)^''^'' which also appear in the Ito derivative of the Z{N) b.c.c. parafermionic 



operator, see Eq. (j49l 



APPENDIX B: NULL VECTORS IN THE PARAFERMIONIC MODULES. 
1. Commutation relations of parafermionic current modes:Zjv sector 



The commutation relations of the mode operators can be deduced from Ea. (|39p by using standard techniques in 
the complex plane. Consider for instance the commutation relation of the modes of ^'^(z) in the Zn singlet sector 
<3?°. We use the short-hand notation j^*^, ^1/1}$" to indicate such relations. One consider the double integral: 



27r« Jco Jco 



^\Z,)^\Z2)^"{0), 



(Bl) 



where one adds the term (zi — ^2)^^^ ^ to make the integrand single- valued. By using the Eqs. (P5|) - ([55|) and expanding 
in series the term (zi — Z2)^/^^^, the Cauchy theorem gives the following relations for {^'^, ^'^}<i>'': 



^{2-N)/N ^(3-JV)/7V-|-n^l/Af+m + ^(3-W)/Ar+m^l 



/JVH 



(4:-N)/N+n+m 



(B2) 



The coefficients are defined from the development 



1=0 



The relations ^}$^, for instance, are derived from the double integral: 



2TTi 



Co 



dzi i dz2Z?/^+"zf-^)/^+"(zi-Z2)-(^+^)/^Vl,l(,,)^-l(,2)ci>l(0). 



Co 



(B3) 



(B4) 



The exponent of the (zi — Z2) term has been chosen on the basis of the OPE In particular it allows the modes 

of the stress energy, appearing at the second order in ([29]) . to enter in the relation One obtains: 



'^D-{N-2)/N ^l/N+n~l^-l/N+m+l ^-l/N+m^l/N+n 



$±1 = 



1 

2 ViV 



— + n - 1 ) (5„+„i,o + 



A^ + 2. 

J 

N 



From the above relations, and using the primary condition p6p . one can readily obtain the conformal dimension Ai 
of a primary field Setting n = to = in Eq.dBS]), one obtains Af = {N - 2)/{N{N + 2)). 

For completeness, we give below in a compact form the the commutation relations j^*^, 5'^}$'' and {^P^, ^'^^}<I>'^ 
which determine the space of the representations : 



(2-N)/N 



(2+N)/N 



1=0 



<j)9 



1.2|g| 



.2|<?| 



TV -1-2 



2^1^ N 



$9. 



(B7) 
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FIG. 5: Representation module of <1>*^ and the action of the parafermionic modes giving the operator e and its second level 
descendants 



The integers s{q), t{q) and u{±l,q) shift the indeces of the parafermionic modes: 

2 

TV ^ 
2 





- "q+2 






- "q+2 


~ K+1 H 


{k,q) 


q 


- 5U - 



1 



N 

N + 2 
N ' 



(B8) 



where (5^_|_2 is defined in Eq. (|37p . Finally the algebra is then completed by the commutators between the parafermions 
modes A^^ and L„: 



^ i„-L„A±i _^ = f(l - A* )n + TO - <5^ , J ^ ^$9.. 



(B9) 



Analogously to what we have seen above for the special case of the field, the formula ^ giving the dimension 
of the primary fields is obtained by taking n = m = (ri = 1, to = 0) in the Ea. (|B7p for g 7^ {q = 0). 

2. module. 



In the following we consider in more detail the module of the representation . The structure of the corresponding 
module is shown in Fig.© for the case N = 5. 

We have seen above that the dimension of the operator are given by the commutation relations Eq. (|B7p without 
imposing any degeneracy in the module. This in general is not the case, as it can be seen in the study of the 
representation of the minimal models or of other parafermionic algebras. However, in order to characterize completely 
the representation one needs to fix the eigenvalue h±i of the zero mode defined in Eq. ([55)) . We show this is easily 
obtained by imposing a degeneracy at level 1/A^ where the first descendent states in this doublet are the two singlet 

Forming the linear combination 



states A_^^j^^~^^ and A^_^^j^ 



=aA\^<i>-^ + bAz\<^\ 



(BIO) 



we wish to make it into a primary operator, i.e., to ensure that it is annihilated upon action by positive index mode 
operators. In this case it will be sufficient to verify that 



A^X°i=0 and A-\x°_±^0. 

^ N N "TlV iV 



(Bll) 
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The degeneracy conditions ^^/V^" i/w ~ ^'"^ satisfied if 



= (B12) 



where the matrix element fj,k,k' is defined by 

. .T>1 — /I „ 



Using Eq.(125D-Eq.(IlSl), one finds: 



2 

^i^i = Xl'^hi, ^1 _i = — , (B14) 



where the structure constant Aj' is given by Eq. (l32|) . and the zero mode eigenvalue hi, defined in Eq. ([38]) . is fixed 
to: 

Once the conformal dimension and the zero eigenvalue Eq. ljBlsp have been determined, the operator x^^±. — 

N 

A^_j_^^^ — A^^i $^ is a primary operator and it is put to zero, x^j_ — Oj to make the representation irreducible. 

N N N 

After imposing the degeneracy at level 1/A^, only one singlet state remains: 

e(^) = (B16) 



with dimension 

2 



N + 2 



(B17) 



3. Independent states at the second level of the e*^^' operator. 

Using systematically the commutation relations Ea. (jB6P -Ea. (jB7p . one can show that at the second level of the e*-^' 
operator there are only two independent states. All the operators at the second level can indeed be obtained by a 
linear combination of the following two states: 

(Alr,A-,\+A-\ A\ ."jeW; and (aI^.A-,^+A-\ ,A\)eW. (B18) 

The action of the above parafermionic modes in the module is shown in Fig.©. 

Now, in order to make the connection with SLE, we express the Virasoro operators L'^i and L_2 in terms of 
combination of parafermionic modes. Using the results Ea. (|B14[ )and Eq.ipB]) in the relations (|B6|) - (|B7p . we can 
express the L-ie operator as: 

L-is^""^ = ^ ^ + A\i/^-i'f-') ■ (B19) 

The above relation, expressed in terms of the coset SU{2)k/U{l), corresponds to the well known Knizhnik- 
Zamolodchikov equations. Starting from Eq. (|B19[) and keeping into account Ea. (|B14p and Eq. ([55|) and other similar 
relations coming from the structure of degeneracies of the module, the commutation relations (|B6P - (jB7p give the 
following relations: 

^-''^"^ = + ^'""^ + ^ {a\/n-2^-' + A-_\,^_,^') (B20) 

L-2e^''' - ^ {az\,,^,A\,,_, + AZI,,,,AI,,^,) + (a}/..^,*- + Az\,,^,<i>^) (B21) 

From the above relation, the Eg. ([55]) is easily obtained. Note that in the case of iV = 3 (the module is not 
present in the N = 2 theory) , the parafermionic algebra provides an additional relation between states at the second 
level of e^^-*. Indeed using the relations {5'^^'^}(f>^^ and the fact that v]/iv|/i — > vp^^, valid for iV = 3, one obtains: 
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Using Eq. (|B22| and Eq. (|BT9l) in Eq. (lB2T1) and in Ea. (lB22l) . one gets: 

i-2 - e(^=^' = (B23) 

which is the relation (|15p obtained by the study of the minimal model M5 (remember that e(^=3) corresponds to (j)i,2 
of the minimal model Kac table). 



4. Commutation relations of parafermionic current modes:i?— sector 

As mentioned in the Section (jlVCp . the reader can refer to and [1^ for a complete discussion about the 
construction of the disorder sector modules and of the derivation of the commutation relations l^*^, ^''^ }Ra- In the 
calculations of degeneracy we have used two types of commutation relations: the first one is between the modes of 
two ^'^ chiral fields and the second one is between the 4"^ and chiral fields, with k = 2, 3, ... , '^^y^. Using the 
expansions Ea. (|43)) and Ea. (|44)) . the {'^'', ^'^} relations have the following form: 



Ra — 



nl-ik/N \ k,l Ak+1 r . n-i+ik/N \k,-l , ixm^fe+l ii 



^k+l^(n+m)/2^ 



Ra' 



(B24) 



where the coefficients D^^ are defined from the expansion: 



(B25) 



and the matrix U'°i?Q — Ra-k changes the R indices. For the {^'^,5'"'^} relations, which from the Eg. ([43]) establish 
the connection with the Virasoro generators, one gets: 



2/N 



^(n-;)/2^(m+;)/2 



^(m-;)/2^(n+0/2 



Ra 



N ■ 



N 



^{n+m)/2 



K(n)Sn 



+m,0 



Ui?a' 



(B26) 



where 



n 



16N 



Finally we will use the commutations {T, 'i'}Ra- 



[Ln,Al^,]= (n+l)A 



m 
2" 



A) 



+m/2- 



(B27) 



(B28) 



The level structure of the modules of disorder operators is relatively simple. There are only integer and half-integer 
levels, and there exists zero modes for all the operators {5'''} acting on the A^— uplet of disorder operators. From the 
expansion (|44)) there are [N/2\ zero modes Aq (with fc = 1, 2, . . . , -^^y^), associated with the parafermion 5''^ which 
acts between the N summit of the module: 



AnRn — hi, U^^Rn 



(B29) 



This defines the eigenvalues hk- We recall that (JRa — Ra-i- The eigenvalues hf, characterize, together with 
the conformal dimension, each representation Ra- Like the case of the Zf^ sector, the formula for the conformal 
dimension together with the zero modes hk are easily derived by solving the system of equations obtained by 
setting n = m = in the Eq. (|B24[ ) and Ea. (|B26p . Here we focus our attention on the module of the operator i?!"' 
with dimension: 



^0 ^ 



N^ + N-2 
16(iV + 2) ■ 



(B30) 
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Rl R2 R3 R4 R5 



1/2 




3/2 



FIG. 6: Representation module of Ra fields for N — 5. The arrows show the action of the parafermionic modes. The zero mode 
action between the summit of the module are illustrated by the blue arrows. 



5. Independent states at the second level of the i?'"' operator. 

We have analyzed more in detail the structure of the module of the _R(°) disorder operator. Putting the values 
n = 0,TO = — 1 in the Ea. (jB24[ ) and Eq. (jB26p . and using Eq. (|B30p together with the value of the zero mode t/i 



2/1 



^-2(N-l)/N^ (B31) 



caractherizing the module of , it can be shown that 

^-1/24°^ -0 k^l,2,..[N/2\ (B32) 

for each k. This simplifies greatly the computation of the relations at the second level of this operator, in which we 
are interested. 

At the first level, by setting n = — 1, m = — 1 in the Eq. (|B26p and using Ea. (|B32p . one gets: 

- 2-^-^/^VnAI,R^I, (B33) 
Combining the above relation and the Eq. (|B28p with n = —1 and m = —2, we find: 

^ 2-i-2/w^^AL2i?i"^i - 2-2-4/A'ivALiAlJi?i"). (B34) 
V N 



Finally, to express the operator L-2 we use Eq. (|B26P with n = —2 and m = —2, obtaining: 

L_,i?i°) = 2-'-y^^Al,R^:i - ^-'""^^AUAZIR^. (B35) 



The Eq. ([M)l is then directly derived from the above relations. 

As in the case of the e^^^'^) operator, we have verified that, in this approach, we find the known results for the (/)2.i 
operator of the A/5 model (remember that for = 3, i?^"^ = 02, i)- Indeed, using the fact that = 'I'^^ for A^ = 3, 
we can extract from the Eq. (|B24p with n ~ —2,m = —2 the following relation: 

A\A-_\R^^^=2y'^A\,R^^\ (B36) 
Using Eq. ([B36l) in Eq.(jB34| and in Eq. ([B35l) . we obtain: 



- ) = (B37) 



in agreement with Eq. (|16p . 
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APPENDIX C: BOUNDARY STATES. 



One of the main results of the Boundary conformal field theory [1^, |l3[ is the bijection between the possible 
conformally invariant boundary conditions and the bulk primary operators. In particular, the allowed boundary 
states are expressed as linear combinations of bulk primary operators. The coefficients of such expansion are directly 
related to the entries of modular transformation matrix S. Reminiscent of the coset construction of the parafermionic 
theories, = SU{2)]\f /U{1), it is convenient to classify the bulk operators with the notation |j, m >, with j = 0, 1, • • • 
and m — — \_N/2\ , \_N/2\ + 1, • • • , \_N/2\ . The boundary states \l,m> are then defined by the formula: 



,m' > 



,77 



s 



0,0 



(CI) 



The three states Potts model represents a typical example where these results apply. In order to fix the notation 
and to facilitate the comparison with the other Ziq theories, we reproduce below the construction of the boundary 
states for this model. We use the same conventions as in [l^. 

The table of the primary fields of the Z-^ theory with the associated characters is: 



Field 


A 


|Z, m > 







|0,0> 


^1 


2/3 


|0,1> 


^-1 


2/3 


|0,-1> 


$1 


1/15 


ll,l > 


£(Ar=3) 


2/5 


|1,0> 


$-1 


1/15 


ll,-l > 



In the basis (|0, >, |1, >, |0, 1 >, |1, 1 >, |0, —1 >, |1, — 1 >) the modular transformation matrix 5 = 5; 
the form: 



I' .m' 



1 



UJS] OJ S] 
9 /' /' 

w s; LOS] 



takes 



(C2) 



where lo = expi7r/3 and 



_ 1 / sin7r/5 sin27r/5 
\/5 I sin27r/5 — sin7r/5 



(C3) 



In order to identify their lattice realization, it is important to study the behavior of the boundary states under a Z3 
rotation. Using the formula (jCl|) . one can observe that under a Z3 rotation the boundary states transform as: 



|0,0> 



10,1 > ^ 10,-1 > 



1,1 > 



1,-1 > 



|0,0 > 
11,0 > 



(C4) 



Identifying the boundary state |0,0 > with the boundary conditions in which all the spin on the real axis take the 
value A, the above transformations suggests the identification of the state |0, 1 > (|0, — 1 >) to the situation in which 
all the spins on the real axis take the value B (C). In consistence with the above transformations, the state |1,0 > 
has been shown to correspond to the state in which the spins can take with equal probability the values B or C . From 
this identification one concludes that the b.c.c. operator ip^ generates the boundary conditions A\B + C discussed 
before. 
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The results shown here for generalize directly to the Zjv theories. The modular transformation properties of 
these theories have been analyzed in [soj . For sake of simplicity we consider here the case Z5. In general, the number 
of irreducible representation of the Virsasoro algebra is infinite for iV > 3. However, there is a system of principal 
Virasoro primary fields which appear at the first levels in each module of The table of the principal primary fields 
of the Z5 theory is the following: 



Field 


A 


Z, m > 


$0 





|0,0 > 


^1 


4/5 


|0,1> 


$-1 


4/5 


|o,-i> 


$2 


6/5 


|0,2> 


^-2 


6/5 


|0,-2> 


$1 


3/35 


|1,1> 


g(Af=5) 


2/7 


|1,0 > 


$-1 


3/35 


|1,-1> 




17/35 


|1,2> 




17/35 


|l,-2> 


$2 


2/35 


|2,2> 




23/35 


|2,1> 




6/7 


|2,0> 




23/35 


|2,-1> 


$2 


2/35 


|2,-2 > 



In the basis (|0,0 >,|1,0 >,|2,0 >,|0, 1 >,|1,1 >, |2, 1 >,|0,2 >,|1,2 >,|2,2 >,|0,-1 >,|1,-1 >,|2,-1 >,|0,-2 > 



,|1,— 2>,|2,— 2>), the modular transformation matrix S = takes the form: 



cm 

Orr, — 



V5 



/ 


*i 




4 






\ 




*i 


wsf 




UJ Si 


UJ 's\ 






*i 




! A J' 

~LJ Si 


UJ ^s\ 


UJ Si 






s'' 
*i 


UJ Si 


UJ ^s\ 


I' 

UJS\ ^ 


UJ s\ 




\ 


s'' 


, i-2o'' 
Ul Si 






^'4 


) 



(C5) 



where uj = expin/b and 



sin7r/7 sin27r/7 sin27r/7 \ 
--y= I sin27r/7 sin47r/7 sin67r/7 j 
sin37r/7 sin67r/7 sin47r/7 / 



(C6) 
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The corresponding boundary states |0,0 > and |1,0 > transforms, under a Z5 transformation, in the following way: 



|0,0 > ^ |0,1 > ^ |0,2 > ^ |0,-2 > ^ |0,-1 > 

|1,0 > ^ |1,1 > -> |1,2 > -> |l,-2 > ^ |1,-1 >. (C7) 

It is clear that the case Z5, and Zn in general, represents, from the point of view of the properties of the boundary 
states, a direct generalization of the case Z3. It is then natural to identify the boundary state |1, > with the state in 
which the spins can take with equal probability the values B, C, D, E. Consequently, the b.c.c operator is expected 
to generates the boundary conditions A\B + C + D + E discussed in section IVl 
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